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Abstract 

Homoclinic bifurcation of 4-dimensional symplectic mappings is asymptotically 
studied. We construct the 2-dimensional stable and unstable manifolds near the sub- 
manifolds which experience exponentially small splitting, and successfully obtain ex- 
ponentially small oscillating terms in the 2-dimensional manifolds. 
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Homoclinic (or heteroclinic) bifurcation plays an important role to cause chaotic motion 
near hyperbolic fixed points. Melnikov's method (or Melnikov's integral), which is based on 
the regular perturbation method, is a powerful tool to detect transversal homoclinic points. 0* 
On the other hand, it is known that in the rapidly forced systems and the standard map, 
and so on, splitting of separatrices is exponentially small. Consequently it is impossible to 
capture exponentially small oscillation straightforwardly with Melnikov's method. 

In the last decade the difficulty has been overcome by using the method called asymp- 
totic expansions beyond all orders. This method enables one to capture exponentially small 
oscillation of the stable and unstable manifolds in 2-dimensional symplectic mappings which 
are perturbed from linear mappings, fl'll'S) and to construct functional approximations of 
the manifolds. i> i> 0> 

For 4-dimensional symplectic mappings, Gelfreigh and Shatarov treated a coupled stan- 
dard mapping and computed the crossing angle between the 2-dimensional stable and unsta- 
ble manifolds, i* We obtained the functional approximations of particular sub-manifolds of 
the stable and unstable manifolds in a 4-dimensional double-well symplectic mapping with a 
weak coupling. S In this letter, we construct the functional approximations of 2-dimensional 
stable and unstable manifolds in 4-dimensional double-well symplectic mappings with more 
general weak couplings. 

We start with the 4-dimensional symplectic map (qj,Pj) (q'j,Pj) 

= p j + e ( qj -2p^-e^ 
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where j, k G {1, 2}, j k, < e 1 is a perturbation parameter. To simplify computation 

we take the coupling potential J as J = —q\q\ H — <?i<?2-^ Note that the map (1) possesses 

2 3 

the substitutional invariance between j and k. 

The map (1) possesses a hyperbolic x hyperbolic fixed point at the origin (0, 0, 0, 0) 
and there exist 2-dimensional stable and unstable manifolds near the origin. In Ref. [|), 
we applied the method to the particular 1-dimensional stable and unstable sub-manifolds 
defined as q\ = q2 and successfully obtained the functional approximations of the sub- 
manifolds when 6 = 0. In this letter we study the neighborhood of the sub- manifolds as a 
perturbation problem. 

One can rewrite the map (1) in the second order difference equation 

A 2 eqj (t) = gj (t) - 2 gj (tf - e\aq 3 {t)q k (tf + bq^tf q k {t)% (2) 



t In Ref. ^) we show that if J is a polynomial of qj 's of degree six, it is possible to capture exponentially 
small oscillation. Hence this coupling term is considerably general. 
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where A 2 e q{t) = {q(t + e) — 2q(t) + q(t — e)}/e 2 . The difference equation (2) is equivalent to 
the infinite order differential equation 

d ~§ = qj ~ 2q i ~ t2 ^ ql + b( t ql) ~ 2 § fol^' (3) 

which is called the outer equation. The boundary condition with which we solve the equation 
(3) is qf(t; e) — > +0 (t — > ±oo), where the superscript ± stands for the stable and unstable 
manifolds, respectively.! 

We first construct the solution to the ODE (3) by using the regular perturbation method 
with respect to e 2 , i.e., expanding qj(t, e) = Sj(t)+e 2 q'j(t) + 0(e 4: ). Substituting the expansion 
into the ODE (3), one obtains the unperturbed equation as 

^fl - g. _ 2S 3 

and the perturbed equation as 

^ = (l-6* i (t) a )fc+F i (t), (4) 

where Fj(t) is a polynomial of Sj(t),Sk{t) and their derivatives. The boundary condition is 
rewritten as Sj(t),q~j(t) —> +0 (t — > ±oo). 

The unperturbed solution is obtained as Sj(t,tj) = sech(t + tj), where t/s are integration 
constants. Because the system is autonomous, it has translational invariance for t. Hence 
the essential parameter is the difference between t% and t 2 . We put tj = (— where 
a is an arbitrary real constant. To simplify notations we put s(t) = Sj(t,0) = sech(t). 
Taking the integration constant a, which was vanished in Ref. as another perturbation 
parameter, we successfully obtain exponentially small oscillating terms in the neighborhood 
of the sub-manifolds. 



The fundamental system of the solution to the homogeneous equation of the ODE (4) is 
given as Sj(t,tj) and gj(t,tj) = g(t + tj), where g(t) = s(t) f * = ^sech(t) — -cosh(t) — 

•/SIT'J — — 
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-t sinh(t)sech(t) 2 . Note that s(t) and g{t) are odd and even functions of t respectively, and 
Wronskian is equal to 1. The general stable and unstable solutions to the inhomogeneous 
equation (4) is obtained as 

Qf(t,tj) = cfs.it) - s 3 (t) f Fj{t) gj {t)dt + gj {t) f FjWsjftdt, (5) 

respectively, where c^'s are integration constants. 



* The same notations will be used without further comment. 
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When a = 0, which is equivalent to tj = 0,j = 1,2, the solution stands for the 1- 
dimensional sub-manifolds of the 2-dimensional stable and unstable manifolds, and the stable 
and unstable sub-manifolds coincide, i.e., q^(t) = qj{— t)M One can therefore restrict the 
solution to even fuctions of t. Putting qf(t, 0) = q(t), one obtains 



~/ x b + 2 , . . o 6a + 86 + 7 , . . t . , , . , , . 9 
g(t) = — — sech(t)' 3 — sech(t) + — sinh(£)sech(t)^. 

Next we construct solutions when tj ^ 0. Let us suppose \tj\ <C 1, corresponding to the 
neighborhood of the sub-maifolds. We expand qj(t,tj) in the power of tj around q(t) as 



qf(t,t J ) = q(t) + J2 t 



dqf(t,tj 



+ 0{t) 

t,=o 



and cf = c° + i,c± + 0(tj). The condition qf(t,tj) -► g(t) (i,- -> 0) implies c° = 0. The 
terms in the bracket of the right hand side of the equation (6) are independent of j because 
of the substitutional invariance between j and k of the map (1). Hence we put them q^(t), 
i.e., qf(t,tj) = q(t) + tjq ± (t). A straightforward computation gives 



q ± (t) = cfs(t) =F 1Q5 + g(t) 



+-6sinh(t)sech(t) 4 J~ — -sinh(t)sech(t) 2 log (cosh(t)) 

8a -66 . , , . , , , 2 28a + 246 . , , , 

H smhmsecnm smhm. 

30 w w 1Q5 w 

Thus the solution to the outer equation (3) is given as 

qf(t,tj) = s(t) + e 2 q(t) +tj{s(t) + e 2 q ± (t)} (7) 

up to e 4 and t 2 . 

The outer solution (7) possesses singular points at t = nij2 + nni, n £ Z. One must 
therefore analyze the neighborhood of the points as the inner problem. Introducing the new 
variables as t — ni/2 = ez, Tj = tj/e and <£>j(z,Tj) = eqj(t,tj), the difference equation (2) is 
rewritten as A 2 <Pj = -2<P 3 j -b<P 2 <I>l + e 2 (<I> j -a<P j <P 2 k ), where A 2 <Pj = <P j (z + l,r j )-2<P j (z,T j ) + 
$j{z — l,Tj). In this letter we restrict our attention to the leading term of <Pj(z,Tj), hence 
truncating the terms with respect to e 2 , we obtain the inner equation as 

A 2 $j = -2<P 3 - - b$ 2 M. (8) 
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Note that when b 7^ 0, the second term, which is generated from the coupling term, contribute 
the inner equation. 

When a = 0, putting X(z) = @j(z, 0), the inner equation corresponing to s(t) + e 2 q(t) is 
obtained as 

A 2 X(z) = -2X(zf - bX(z) 5 , (9) 

which is called the first inner equation. The matching condition is X(z) — > e(s(i) + e 2 q(t)) 
when \z\ — > oo. On the other hand, when a^O, expanding ^(t, Tj) = X(z) + TjX'j(z), we 
obtain 

A 2 X;(z) = -QXtzfX'^z) - bX(zf{2X' 3 {z) - 3X^)}. (10) 

Summing up both sides of the equations (10), one obtains the second inner equation 

A 2 X(z) = -QX(z) 2 X(z) + bX(z) 4 X(z), (11) 

where X(z) = X[(z) + X^(z)^ The matching condition is X{z) — > e 2 (s(t) + e 2 g ± (t)) when 
\z\ — > oo. 

The inner equations (9) and (11) possess formal solutions when \z\ tends to oo. Although 
the formal solutions do not converge, the formal Borel transformations of these formal solu- 
tions converge, i.e., these formal solutions are Borel summable. By using Borel transforma- 
tion, Stokes phenomena, which are abrupt jump of asymptotic expansions, around t = m/2, 
can be captured, and the analytic continuation from the left half of t-plane to the right one 
can be performed. The details of the computation will be reported elsewhere. 

Finally the stable and unstable solutions which are valid both the left and right half 
plane are obtained as 

#(t,t j ) = tf(t)+t j qt(t), (12) 

where 

Qo{t) = s(t) + e 2 ~q(t) + S(t, e)M(t, e)B(t, e) 

qf(t) = s(t) + e 2 g±(t) + S(t, e){M(t, ej^t, e) + M(t, e)B 2 (t, e)} (13) 
M(t,e) = ---e~- 

M(t,e) = -- 7 re-r 

§ The subtraction of both sides of the equations (10) gives ±he linearized equation of the equation (9). 
In this computation the linearized equation has no information. E3* 
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„/ N 3 2irt . 2nt 

B(t, e) = -TTs(t) cos h g{t) sm 

4 e e 

/ \ 3 .. . . 27rt . 27rt 

-cMt, e) = -TTs{t) cos h gm sm 

4 e e 

fw \ 3 . . . . 27rt 27rt 

±32(^7 e) = "jTrsf t) sm gm cos 

I e e 

(t < 0) 

e) = <( 1 (t = 0) (for small e). 

2 (t > 0) 

The Stokes constants c and c, which depend on b, are numerically obtained as c = —763 and 
c = 1146 when b = 1.0. For other values of b the Stokes constants also converge. 1 

By using the approximate solution (12), the splitting angles between the stable and 
unstable manifolds can be computed. One can show that the angle along the sub-manifolds 
(ti = £2) is exponentially smallB* and the angle orthogonal to the sub-manifolds are 0(e 2 ). 

The expression (13) implies that not only the sub-manifolds but also the 2-dimensional 
stable and unstable manifolds themselves possess exponentially small oscillating terms. In 
the future work, we extend this computation to 2n-dimensional symplectic mappings with 
weak coupling terms. 
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